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In a hybrid optomechanical setup consisting of a two-level atom in a cavity with a pendular 
end-mirror, the interplay between the light field’s radiation pressure on the mirror and the dipole 
interaction with the atom can lead to an effect, which manifests itself in the suppression of Rabi 
oscillations of the atomic population. This effect is present when the system is in the single-photon 
strong coupling regime and has an analogy in the photon blockade of optomechanics. 

PACS numbers: 42.50.Wk, 42.50.Pq, 37.30.-|-i 


I. INTRODUCTION 

In the Jaynes-Cummings model, a single atomic dipole 
interacting with a single photon of an electromagnetic 
field mode in a cavity periodically exchanges the excita¬ 
tion between the electronic and photonic degree of free¬ 
dom. In this letter we ask ourselves: How does this be¬ 
havior change, when the cavity’s mirrors are not fixed, 
but can move in time as an additional quantum dynam¬ 
ical degree of freedom? How does the quantumness of 
the mirror motion manifest itself in the dynamics? And, 
what happens if several photons are present in the cavity 
field? 

The recent experimental success in cavity optomechan¬ 
ics , notably the demonstration of laser cooling of the 
cavity’s mirror motion towards the ground state of the 
confining harmonic potential, pioneers near-future exper¬ 
iments with extended, hybrid systems. An obvious ex¬ 
tension of the paradigmatic optomechanical system, i.e 
an optical cavity coupled to a mechanical element by ra¬ 
diation forces whose motion changes the cavity’s bound¬ 
aries, can be achieved by adding a single two-level sys¬ 
tem which couples to the cavity light field by dipole in¬ 
teraction, thereby combining cavity optomechanics with 
cavity quantum electrodynamics. Besides conventional 
quantum electrodynamics experiments with macroscopic 
cavities investigations of such systems are of particular 
importance for solid-state systemSjSuch as two-level sys¬ 
tems in optomechanical crystals (a-Q- In circuit elec¬ 
tromechanical setups the coupli ng to two-level systems 
has already been demonstrated [§-[lll|. 

Especially solid-state systems (leaving aside the two- 
level system in contemporary realizations) are candidates 
for reaching the so-called single-photon strong coupling 
regime [l^, where the radiation force of a single photon 
is strong enough to displace the harmonically oscillating 
mechanical element by more than the extension of its 
ground state wave packet. This regime is entered, when 
the mechanical frequency becomes much smaller than the 
optomechanical coupling. Then, the non-linearity of the 


* ralf.betzholz@physik.uni-saarland.de 


radiation force becomes clearly apparent and manifests 
itself in phenomena like the photon blockade effect [l3l| . 
which is caused by an effective photon-photon interaction 
mediated by the mechanical degree of freedom. 

When adding a single atom to the optomechanical cav¬ 
ity, the resulting hybrid system does not only reflect the 
radiation forces of the cavity-mirror interaction, but also 
includes the dipole interaction between a near-resonant 
atomic transition and the cavity field, drastically altering 
the dynamics as a whole. The richness in the dynamics 
of such a highly non-linear system has been started to 
be explored in several general investigations of its co¬ 
herent motion [l3 - [l^ and can be exploited to prepare 
non-classical states, similar as in nonlinear optical se¬ 
tups [ 13 , [ 11 . In this letter we focus on the possibility to 
dynamically suppress the Rabi oscillation of the atomic 
population due to the mirror motion. The idea behind 
can best be illustrated having a classical picture in mind, 
i.e. when the mirror is moving parametrically along a 
given trajectory. Then, an initially tuned atom-cavity 
interaction becomes disturbed by the mirror motion due 
to the change of the cavity’s resonance frequency being 
a function of the cavity length. If the mirror elongation 
is large enough, the modified atom-cavity detuning can 
lead to reduced Rabi oscillations of the quantum electro¬ 
dynamic subsystem. We will develop this simple classical 
picture beyond the trivial parametric motion towards a 
fully dynamical one in the quantum regime. The Wigner 
function will be used to represent the dynamics of the 
mechanical degree of freedom, which can exhibit strong 
non-classical characteristics deep in the non-linear pa¬ 
rameter regime. Even when dissipation and variation of 
the initial state is included, characteristic features of the 
suppression can be distinguished. 


II. THE HYBRID OPTOMECHANICAL SETUP 

We consider an optomechanical resonator with a two- 
level atom placed at a fixed position within the resonator. 
The transition frequency between the ground state \g) 
and the excited state |e) is assumed to be resonant with 
the frequency ca of the cavity when its end mirror is in 
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its equilibrium position. The atomic dipole couples via 
Jaynes-Cummings interaction [l^ to the electromagnetic 
field of the cavity while the cavity interacts with the me¬ 
chanical oscillator of frequency v and mass M by radi¬ 
ation pressure. We denote the dipole-cavity interaction 
strength and the optomechanical coupling strength by g 
and x^ respectively. A schematic picture of this hybrid 
quantum system is shown in Fig. [TJ The total Hamilto- 


with the definition /3 = and displacement operators 
D{a) = exp[a6l — a*b]. The states \n) and |m)mec stand 
for the cavity and mechanical Fock states, respectively. 
The propagation of the mechanical oscillator, initially in 
prepared in its ground state, is given by the coherent 
state [2l| 

|»7(i))mec = |n/3(l - e“*''‘))inec ■ (6) 



FIG. 1. Hybrid optomechanical setup consisting of a single¬ 
mode cavity with a moving end-mirror confined in a harmonic 
potential of frequency v and a resonant two-level atom at a 
fixed position. The cavity couples to the mechanical oscillator 
by radiation pressure with a coupling strength x S'lid to the 
two-level atom via dipole interaction of strength g. 

nian of the hybrid systems reads 

= iJom + iJtls + iJjC (1) 

where the operator iJom labels the standard Hamiltonian 
of optomechanics 

iJom = fuj-ia^a + hvb% — (b -|- . (2) 

Here, a and represent annihilation and creation oper¬ 
ators of the cavity mode, and b and b^ for the mechanical 
oscillator. Its position operator relative to the equilib¬ 
rium is X = ^(&-l-6^) where the harmonic oscillator length 
scale ^ = yJ'h/2Mv denotes the extension of the ground 
state wave packet. The Hamiltonians of the two-level 
system and the Jaynes-Cummings interaction are given 

by 

i/tis = /ia;|e)(e|, (3) 

^JC = ^5 [a'^\9){e\ + a\e){g\] , (4) 

respectively. 


III. DYNAMICS OF THE OPTOMECHANICAL 
SUBSYSTEM 


Figure (St-) shows the trajectories in Wigner phase space 
for the values n = 1 and n = 2 of the cavity photon 
number. The radiation pressure pushes the mirror out 
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FIG. 2. a) Gircular trajectories of the mechanical oscillator’s 
Wigner function for cavity photon numbers n = 1 (solid) 
and n = 2 (dashed) if the mirror was initially prepared in 
its ground state. The radii are proportional to /3 while for 
n = 0 the Wigner function rests centered around the origin, 
b) Trajectory of the oscillator wave packet in the harmonic po¬ 
tential modified by the radiation pressure. The (a;)max mark 
the maximal elongations. 

of its equilibrium position, resulting in a modified cavity 
length and thereby a modified frequency. The maximum 
displacement from the equilibrium position 

= 2^ Re[n/3(1 - = 4n/3^ (7) 

is reached at time t = tt jv. In the single-photon strong 
coupling regime this displacement amounts to a signifi¬ 
cant change in the cavity frequency during the dynamics. 


Before we delve into the motion of the full hybrid sys¬ 
tem we set the stage by recalling the dynamics of the op¬ 
tomechanical system in absence of the atom. The eigen¬ 
states of Horn are given by [2l| 

(5) 


IV. COHERENT DYNAMICS OF THE HYBRID 
SYSTEM 

For the coherent dynamics of the hybrid system in¬ 
cluding cavity, atom and mechanical oscillator, we focus 
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on the initial state = 0)) = |e)|0)|0)i„ec- Exper¬ 
imentally such a state can be prepared with the help 
of ground-state cooling of the mechanical oscillator and 
standard optical pumping for the atomic state. Ground- 
state cooling for an optomechanical crystal was already 
demonstrated and is expected to be realizable in sim¬ 
ilar hybrid systems in near-future experiments. 

The position of the mirror determines the frequency of 
the cavity, which comes out most clearly when rewriting 
the Hamiltonian ([T]) in the form 


H = h 


X 

U! — —X 


i^a + hi/b^b + Htis + Hjc ■ (8) 


The mechanical oscillator in its maximally elongated 
mean position Eq. o results in an effective detuning 
Sn = 4n/3x between the cavity and the atom. From this 
consideration follows that for a distinct effect of the mir¬ 
ror’s motion on the Rabi oscillation one has to be in the 
strong-coupling regime with /3 ^ 1. In this regime we 
cannot resort to the a ppr oximate solution of single po- 
lariton optomechanics j22l |. 


A. Limiting cases 

Nevertheless, to provide first insight into the dynamics 
we start our investigation with two limiting cases, namely 
g <C X, and /3 -C 1, that allow for an approximate 
analytical treatment. 


1. Slow Rabi oscillations 


In the case of small Jaynes-Cummings coupling, i.e. 
g focus on the parameter region /3 Ri 1 

for strong optomechanical coupling and restrict the fol¬ 
lowing discussion to the subspace of only one excitation 
in the quantum electrodynamical subsystem, remarking 
that the generalization is straightforward. In the inter¬ 
action picture and after displacing the Hamiltonian o 
with D{—P\g, l){g, 1|) it takes on the form 


H = hg 


e*'’'^‘ZI(/3e*‘^*)|e,0)(g,l| -f H.c. 


(9) 


apart from a constant energy, in the relevant subspace 
spanned by the two states |(?, 1) = |g)|l) and |e, 0) = 
|e)|0). A rotating-wave approximation can be performed 
after expanding the displacement operator in a power 
series and taking only those terms that rotate with the 
smallest occurring frequency, which for /3 Ri 1 is given 
by — l)i^. In the expansion of D{f3explwt]) these 
surviving terms contain a single operator b more than b^. 
Applying this procedure to (0 leads to the approximate 
Hamiltonian 


where the function /(to) can be obtained by evaluating 
the matrix elements niec(TO|iJ(/3)|TO -|- l)niec [13 which 
explicitly results in the expression 

with the generalized Laguerre polynomials Lm\x). 
Within this approximation the Hamiltonian is a Jaynes- 
Cummings-type interaction between the atom-cavity 
subsystem and the mechanical oscillator with an addi¬ 
tional phonon number dependent factor that does not 
alter the eigenstates, i.e. the dressed states. In the orig¬ 
inal picture the Hamiltonian then reads 

^Irwa =hiyb'<b-nx\g,l){9A\ib + b'‘) 

+ hg [f{h%)hD\p)\e, 0)(5,1| + H.c.] (12) 

and its eigenstates are dressed states in which the 
state Ig, l)|TO)i„ec is displaced by U(/3) while the state 
|e, 0 )|to -I- l)mec remains undisplaced. For the initial 
state where the atom is prepared in its excited state 
and the mechanical oscillator in its ground state, i.e. 
\'ip(t = 0)) = |e, 0)|0)mec, the time evolution results in 
the reduced density operators 

AtRWA(t) = [cOS^(Ht) -|-COS^1?sin^(Ht)] |0)niec(0| 

-I- sin^'dsin^(Ht)£)(/3)|l)i„ec(l|T’l(/3) (13) 

and 

PRWA(t) = [cos^(Ht) -|-cos^ i?sin^(Ht)] |e,0)(e,0| 

+ sm‘^'dsm^{nt)\g,l){g,l\ (14) 

for the mechanical oscillator and the atom-cavity subsys¬ 
tem, respectively. Here we defined tani? = — l)v 

and H = 2gf{0). We compared these approximate re¬ 
sults with numerical propagation and found especially 
good agreement for the case /3 = 1. From these out¬ 
comes we conclude: (i) The motion of the mechanical 
oscillator leads to non-classical states being composed of 
two contributions, namely an incoherent superposition of 
the vacuum and a displaced Fock state. They correspond 
to the radiation pressure of zero and one photon which 
are associated with the two trajectories n = 0 and n = 1 
of Fig. [2^) in the uncoupled case and become mixed due 
to the coherent Rabi oscillations between |g, 1) and |e, 0). 
(ii) For (3^1 one finds detuned Rabi oscillations of the 
atom-cavity subsystem, while they are strictly sinusoidal 
for /3 = 1. (Hi) The regime considered here is not suffi¬ 
cient for observing a suppression of Rabi oscillations of 
the atomic population. Nevertheless, it provides a first 
insight into the coupled dynamics of the system and the 
origin of the emergence of non-classicalities in the state 
of the mechanical oscillator. 


2. Small optomechanical coupling 


HnwA = hg [e*(^'-b-V(6t6)6|e,0)(5,l| +H.C.1 (10) 


Another limit which allows for an analytical treatment 
is the regime of small values of (3 where Ref. [13 provides 
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an approximation of the Hamiltonian ([I}. According to 
that treatment, in a frame displaced by D{—P/2) with 
/3 <C 1 and apart from a constant term, the dynamics is 
governed by the Hamiltonian 

Hp =hvh^b + hg [|+)(+| - |-)(-|] 

- ^hi^[b\+){-\+Il.c.] , (15) 

in the subspace spanned by the usual Jaynes-Cummings 
dressed states |±) = [|(?, 1) ± |e, 0)]/-\/2 of the resonant 
atom-cavity subsystem. Eq. (USD is again of Jaynes- 
Cummings form, with the ladder operators |±)(t| play¬ 
ing the role of the atomic rising and lowering opera¬ 
tors and the mechanical oscillator taking over the cav¬ 
ity part. The approximate form (TTSl) is based on a ro¬ 
tating wave approximation which additionally requires 


\2g — i/l 2g + ly. In the case of a resonant interaction, 
2g = ly, the double dressed states 


|0,0) =|-)|0)„,ec 

|±,m) ~{\+)\m 


(for m = 0), 

l)mec i I )|^)mec] 


(16) 

(17) 


are the eigenstates of the Hamiltonian (fTSl) with the cor¬ 
responding eigenfrequencies Eq = —v 12. and 


£ 


(±) = 


w ~ X- T w 
2 2 



V 


(18) 


for m > 1. The initial state = 0)) = |e,0)|0)i„ec can 
be propagated with help of the Hamiltonian (USD and 
leads to 


li’it)) = e ^ 


m—0 


•\/2m! 


- isin yy/m^iytj e"'^*\+)\m - l)i„ec 
-I- cos e“*^*|-|-) — cos ) 

-I- isin ^\/niTT^iyt^ -|- l)mec| 


\m)r 


(19) 


in the displaced frame. More insight reveals the m = 0 
term. 


IV'(i)) 


e 8 


cos( )e-*t‘|+)-e*t‘|-) 


\l3/2)r 


+*sin(^^z.t^e-*t‘|_)Z5(|)|l)^^^| , (20) 


written here in the original frame, showing the domi¬ 
nant behaviour of |'0(t)) for /3 <C 1: The first term de¬ 
scribes the usual Rabi-oscillations of the atom-cavity sys¬ 
tem when ,8—^0, whereas the second term represents an 
oscillating displaced number state of the mechanical os¬ 
cillator, associated with the dressed state |—). Neverthe¬ 
less, the condition /3 -C 1 for the validity of the Hamilto¬ 
nian (fTSl) prevents clear signatures of non-classicality of 
the evolved state in this limit. 

The time evolution of the excited state population 
Pe{t) = Tr[|e)(e|i/;(t))('0(<)|], using the state (I20l) . can be 
further approximated by 


Pe{t) 


1 

2 


1 -I- cos 



cos(i^t) 


( 21 ) 


which constitutes a sinusoidally modulated Rabi oscilla¬ 
tion with a beat frequency of jiv 12 which describes the 
onset of the suppression of Rabi-oscillations for small val¬ 
ues of fi. 


B. Suppression of Rabi oscillations 

We now focus on the strong-coupling regime ^ 1 
where the approximate Hamiltonians are not valid and 
one has to deal with the hybrid Hamiltoniaii, Eq. o, 
that can be mapped onto a driven Rabi model [24| . Since 
such a model does not provide easy-to-handle analytic 
solutions [1^, [ 2 ^ we restrict ourselves in the following on 
the discussion of numerical results. 


1. The case g = v 

We first focus on the case g = v and /3 = 1. We start 
again from the same initial state = 0)) = |e, 0)|0)inec 
and propagate it numerically using the Hamiltonian ©■ 
In Fig. 13] we show the time evolution of the excited state 
population Pe{t) = Tr[|e)(e|'0(<))('0(t)|]. The strictly pe¬ 
riodic Rabi oscillation of the uncoupled case (dashed) un¬ 
dergoes a drastic change. Instead of the sinusoidal time 
evolution the excited state population exhibits a strong 
suppression before it rises up again at a time appreciably 
longer than the Rabi period tt/^. This behavior continues 
quasi-periodically for the considered parameters and can 
be qualitatively understood with the following intuitive 
explanation: As time evolves the atom initially prepared 
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t/r 



FIG. 3. a) Time evolution of the excited state population 
for the initial state |e, 0)|0)mec. The population experiences 
a strong suppression before a first revival at a time larger 
then a Rabi period. The parameters are /3 = 1 and g = v. 
For comparison the unperturbed Rabi oscillations (/3 = 0) 
are shown in the dashed line, b) Inset on a time scale until 
the hrst revival occurs. The circles indicate instances of time 
which will be referred to later. 


in the excited state populates the cavity with a single 
photon leading to a rising radiation pressure inside the 
cavity pushing the mirror significantly out of its equilib¬ 
rium position and driving it along the trajectory n = 1 
in Fig. [2^). Its displacement is accompanied by an effec¬ 
tive dynamical detuning between atom and cavity sup¬ 
pressing the Rabi oscillations. After approximately one 
mechanical oscillator period r = 27r/z/ the atom-cavity 
resonance condition is fulfilled again and the Rabi oscil¬ 
lation continues. This simple picture however does not 
explain all details in Fig. The behavior is not strictly 
periodic in the mechanical period r and Pe{t) does not 
show a full revival. 

To systematically analyze the suppression of the Rabi 
oscillations for different values of (3 using the same ini¬ 
tial state \'ip(t = 0)) = |e, 0)|0)i„ec as before we show in 
Figure 0^) the time evolution of the population of the 
excited state in dependence of /? color-coded in a den¬ 
sity plot for g = v where dark shading indicates high 
values of the population in |e). For the dotted white 
lines in the lower part of Fig. 0^) we used the pertur¬ 
bative results of Sec. IIV A 21 to show the times of the 
n-th maximum of the perturbed Rabi oscillation for the 
borderline parameters chosen here. These times can be 
approximated by the expression r„ Ri n7r[l — mn/3^/2]/5, 
where the play the role of average phonon numbers. 
As P increases this quadratic behaviour ceases and dis¬ 
tinct maxima of the population emerge separated by gaps 
approximately given by the mechanical period r being 



t/r 

b)d = l 





0 1 2 


x/2/3^ 


c) /3 = 2 



- 


- 

t^0.4r 



t = 2r 
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FIG. 4. a) Dependence of the Rabi oscillations on the op¬ 
tomechanical coupling ratio P for g = v. The dotted lines 
show the dependence of the times at which the population 
exhibits maxima in dependence of /3 following from a pertur¬ 
bative treatment of the optomechanical coupling, b) Wigner 
function of the mechanical oscillator for the eight instances 
denoted by the circles in a) along the line /3 = 1. c) Wigner 
function of the mechanical oscillator for the eight instances 
denoted by the circles in a) along the line /3 = 2. Encir¬ 
cled areas denote regions where the Wigner function takes on 
negative values. 


significantly longer than a Rabi period r/2 for the pa¬ 
rameters. For even larger /3 ^ 3 the first revival sta¬ 
bilizes while subsequent revivals of the population show 
interference-like patterns and decline. 

In order to provide deeper insight into this behaviour 
we show in Figs. 0 }d) and c) the Wigner function of the 
mechanical oscillator for two different values of /3 and 
eight different instances of time as marked by the open 
circles in the density plot Fig. 0^) (and in Fig. 0 ]d) for 
/3 = 1). For both cases the mechanical oscillator starts 
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in the vacuum state represented by the Gaussian Wigner 
function of minimal uncertainty. For (3 = 1 the phase 
space distribution begins to follow the trajectory belong¬ 
ing to the radiation pressure of a single photon inside the 
cavity until at t ~ 2r/5 the quasi probability distribution 
splits into two contributions and the non-classicality of 
this state manifests itself in the negative parts in between 
which are denoted by encircled areas. This behavior is 
caused by the coherent swapping between the states |e, 0) 
and 1 5 ,1) of different radiation pressures corresponding 
to the two trajectories n = 0 and n = 1 in Fig. [5^). At 
t Ri 3r/5 the two contribution are clearly distinguishable 
and merge again towards the end of the mechanical pe¬ 
riod at t ~ T hence bringing atom and cavity back into 
resonance. One can observe however that the Wigner 
function is mostly localized around t 6t/5, i.e. after 
more than one mechanical period. We conjecture that 
this delay stems from the Rabi oscillations taking place 
simultaneously with the oscillation of the mechanical ele¬ 
ment. After that, the quasi-periodic behavior starts into 
another round. For (3 = 2 the time evolution of the 
Wigner function is coined by the appearance of strong in¬ 
terferences in phase space for almost all instances of time 
as shown in Fig. [dj;). Only at < = 0 and the first revival 
at t ~ 1.3t it is positive, already at the next revival at 
t Ri 2.6r negativities and a delocalized Wigner function 
can be observed. This spread in phase space and the in¬ 
terferences explain incomplete population of the excited 
state and the fading of the later revivals: no sharp mirror 
position and thereby no defined resonance between atom 
and cavity can be found for later times. We also checked 
if the blockade effect persists for higher photon numbers 
and found indeed very similar characteristics. 


2. The case g = v 12 

We briefly consider the case g = ;//2, whose limit for 
small /3 was analyzed in Sec. IIV A 21 The numerical re¬ 
sults in this case for /? = 1 are presented in Fig. [5] The 
excited state population in subfigure a) alternatingly dis¬ 
plays incomplete and almost full revivals separated by 
roughly At = 1.75t. In contrast to the previous case 
g = V, the Wigner functions of the mechanical oscillator 
shown in Fig. [SId) exhibit much more pronounced neg¬ 
ativities. This phase space structure of the mechanical 
oscillator can be interpreted as an extrapolation of the 
simple superposition of vacuum and Fock state |l)mec, 
Eq. (1^ . being present for /3 —>■ 0, towards (3 = 1. 

Before we go over to the dissipative dynamics we point 
out the connection of our model to the photon blockade 
in optomechanical systems. The photon blockade [l^ re¬ 
sults from the interaction between the cavity field mode 
and the mechanical object due to the radiation pressure 
which lifts the degeneracy of the equidistant spectrum 
of the harmonic oscillators. Descriptively expressed, if 
a resonant photon is scattered into the cavity, a second 
photon is effectively detuned from resonance and there- 



FIG. 5. a) Time evolution of the excited state population for 
the initial state |e,0)|0)mec for (3 = 1 and g = v(2. The un¬ 
perturbed Rabi oscillations (d = 0) are shown in the dashed 
line, b) Wigner function of the mechanical oscillator state at 
the eight instances denoted by open circles in a). The encir¬ 
cled areas indicate regions where the Wigner functions takes 
on negative values, more circles correspond to more negative 
values. 


fore reflected with increased probability. When the enter¬ 
ing photon stems from a second cavity which coherently 
couples to the optomechanical resonator in the form of a 
beam-splitter-like interaction, the analogy to the model 
presented is established by replacing this second cavity 
with a two-level atom. 


V. DISSIPATIVE DYNAMICS 


In this section we include spontaneous decay of the 
atom with rate F, cavity losses at rate k and damping 
of the mechanical oscillator on a time scale 1 /y into ac¬ 
count. In the regime considered in this work the dissipa¬ 
tive dynamics can be described by the master equation 






( 22 ) 


in Born-Markov approximation. The strong optomechan¬ 
ical coupling requires a treatment where the damping 
of the cavity and the mechanical oscillator cannot be 
treated independently ^i\ . [^ . The non-unitary parts 
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of the dynamics are given by 

C-t\sQ =^'D[\g){e\]Q, (23) 

7 rji 

^cavQ = ^'D[a]g + 4'y-^'D[Pa^a]g , (24) 

CmecQ = '^rh'D[b'^ - I3a^a]g+ + l)75[6 - (3a^a]g 

(25) 

where we introduced the super-operator notation 
T>[X]g = 2XgX’^ — X’^Xg— gX^X to describe Lindblad 
terms. The mean thermal phononic occupation at tem¬ 
perature T is given by m = [exp(?ij^/fcBT) — 1]”^. Again 
the same initial state g{t = 0) = |^/’(t = O))('0(t = 0)| was 
chosen and its time evolution was calculated by numerical 
diagonalization of the Liouville operator £ in Eq. . 

When the dissipative dynamics is dominated by the 
mechanical damping, i.e. 7 ^ k, E, the mechanical oscil¬ 
lator approaches its thermal state and thereby averages 
out the Rabi oscillations of the atom-cavity subsystem to¬ 
wards a constant value of the atomic population as shown 
in Fig. | 6 ^) for 7 = Q.\v (solid) and 7 = OAi/ (dashed). 
Nevertheless, the suppressed Rabi oscillation are observ¬ 
able even for such mechanical decay rates which are much 
larger than in typical existing or planned optomechanical 
setups with Q-factors in the order of 10^ — 10® [2^-[3lj. 

For high-Q mechanical oscillators the relevant time 
scale for damping is given by the atomic and cavity decay. 
Such a situation is depicted in Fig. [BJo) for k = F = O.lz/ 
and (3 = 2 (solid) where we also show the undamped case 
(dashed) for comparison. For such high values of /3 where 
the Wigner function exhibits pronounced negativities one 
might suppose that the virtually undamped mechani¬ 
cal oscillator is capable to support a non-classical quasi¬ 
stationary state. However, the snapshots in Fig. | 6 }:) re¬ 
veal that the time evolution disembogues into a non¬ 
trivial quasi-stationary Wigner function where however 
all negativities are averaged out due to the randomness 
of the spontaneous emission and cavity losses. We re¬ 
mark that in an experiment where the emitted photon 
is recorded the negativity of the Wigner function persist 
for such a single quantum trajectory [s^. 


VI. INITIAL STATES 


We eventually analyze the dependence of the modified 
Rabi oscillations on the initial state. Current technology 
provides techniques which allow for the preparation of 
atomic and cavity states to a high degree of fidelity. The 
initial state of the mechanical oscillator is prepared by 
laser cooling. Most relevant initial states to be considered 
are hence displaced thermal states of the form 


/^a,m — '^(^) 


1 

/ _ \b^b~ 

1 m \ 

771-1-1 

\7Tl - 1 - 1 ) 


D^a), (26) 



x/2/3e 


FIG. 6. Time evolution of the excited state population under 
dissipative dynamics for g = v. a) The bad oscillator case 
(r = K = 0) for d = 1- The mean phononic occupation 
is m = 0.25 and the mechanical decay rates are 7 = Q.lv 
(solid) and 7 = 0.4;^ (dashed), b) The good oscillator case 
(7=0) with the atomic and cavity decay rates F = k = O.W 
(solid) and the undamped case (dashed) for {3 = 2. c) Wigner 
function of the mechanical oscillator for the good resonator 
case of b) at times t = r, t = 4r and t = 50r. Encircled areas 
denote negativities of the Wigner function. 


such that g(t = 0) = |e, 0)(e, 0\g,a,ffi- The initial mechan¬ 
ical displacement is characterized by the complex value 
a and the extension in phase space is proportional to the 
mean thermal occupation fh. The variance of the me¬ 
chanical number operator in these states has the value 
Am^ = fh{rh -|- 1) -I- (2m -I- l)|Q!p. We restrict ourself to 
the familiar case oi g = u and (3 = 1. In order to compare 
the behavior of the modified Rabi oscillations with the 
ideal case, i.e. an initial mechanical vacuum, we adopt 
the measure 

^6r/5 

F(a,m) = / At , (27) 

representing the mean quadratic deviation of the atomic 
population belonging to from the ideal 

case integrated over a time window where the 

suppression takes place. In Fig [Tja,) we show a density 
plot of F{a,m) for |a| < 1 and m < 2 in the good res¬ 
onator case with k = T = 0.05iz and coupling strengths 
g = u and (3 = 1. For the complex parameter a we 




























FIG. 7. a) Density plot of the measure F{a, m) for the sup¬ 
pression of the Rabi oscillation in dependence of |a| and fh 
of the initial state (for each modulus |a| an average over 
eight equidistant polar angles was performed). The dashed 
lines indicate the constant values of the phonon number vari¬ 
ance Am^ = 1,2.75,6. b) Shows two exemplary evolutions 
of the excited state population for the initial states denoted 
by open circles in a) with |a| = 0.5, m = 1 (solid line) and 
|a| = 1 , m = 2 (dashed-dotted line) corresponding to the 
values F = 0.15 and F = 0.25, respectively. As references 
the case of the usual initial state, viz. a = 0 , m = 0 , for 
K = r = 0 is shown in the dashed line and the unperturbed 
Rabi oscillations in the dotted line. The parameters are g = v, 
/I = 1, 7 = 0 and re = F = 0.05i^. 


preformed an average over eight equidistant polar angles 
for each modulus |q;|. From the density plot a quite sta¬ 
ble behavior of the suppressed Rabi oscillations can be 
read off, which only marginally diminishes for tempera¬ 
tures corresponding to m < 0.25 and displacements up 
to |q;| = 0.5, but even higher temperatures and displace¬ 
ments are tolerable. 

To clarify the expressiveness of the deployed measure 
we additionally depict the time-evolution of the atomic 
population in Fig[7lD) for the different initial states indi¬ 
cated by open circles in Fig. [7^), i.e. |a| = 0.5, fh = 1 
and a = 1, 771 = 2 which correspond to F = 0.15 and 
F = 0.25, respectively. Strongly displaced thermal states 
hence do not exhibit an appreciable maximum in the re¬ 
gion around t = 6t/ 5, as found in the ideal case (dashed 
line and Fig. |^. The higher the displacement and mean 
thermal occupation the flatter the curves for the time 
evolution of the excited state population become, thereby 


drastically diminishing the visibility of the suppression 
effect. Nevertheless, within a quite large region in the 
parameter space of the initial states, the effect remains 
clearly distinguishable. 


VII. CONCLUSION 

We identified a distinctive feature in the highly non¬ 
linear regime of a hybrid optomechanical system which 
manifests itself in the suppression of Rabi oscillations 
of the atom-cavity subsystem and dynamically produces 
non-classical Wigner functions in the mechanical degree 
of freedom. The presented phenomenon is linked to 
blockade effects of non-linear systems, particularly to the 
photon blockade of optomechanics. We analyzed the dy¬ 
namics of this suppression and gave a qualitative and 
intuitive explanation using the time evolution of the me¬ 
chanical Wigner function. We also pointed out that even 
under the influence of dissipation the main features of 
the effect are still observable. 

For the case of a resonant atom-cavity interaction, the 
effect is clearly observable when 13 > 1 and g ~ v- Since 
it is based on the coherent time evolution, strong cou¬ 
pling is required, meaning re, F <C g and 7 zz. The 
latter requirement is already achieved in the majority of 
optomechanical experiments with high-Q mechanical ele¬ 
ments. Ground state cooling is required such that the ini¬ 
tial state fulfills < 1. We exemplify these requirement 
based on current experiments, such as the experiment of 
Ref. [ll, with a mechanical frequency v = 2 'k ■ 4 GHz, 
mechanical (optical) Q-factor of Q = 10® {Q = 10®) and 
fh K, 1, only lacking the necessary strong optomechani¬ 
cal coupling which should, along with the cavity-dipole 
coupling, also be in the GHz range. In upcoming hybrid- 
systems with diamond-based crystal cavities [a 1 M|, a 
resonant interaction with color centers with linewidths 
in the MHz range could be tailored, basically fulfilling 
the requirements for the observation of the suppression 
effect presented here. 

With this work we provided deeper insight into the 
rather unexplored area of the non-linear dynamics of hy¬ 
brid optomechanical systems being potentially pivotal for 
future quantum technological applications interfacing dif¬ 
ferent quantum degrees of freedom. 
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